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Abstract 

We show that, contrary to the claim made in arXiv:0911.1299, the extended Hofava 
© : gravity model proposed in arX 1 v:0909.3525 does not suffer from a strong coupling prob- 

I/"} . lem. By studying the observational constraints on the model we determine the bounds 



on the scale of the ultraviolet modification for which the proposal yields a phenomeno- 
logically viable, renormalizable and weakly coupled model of quantum gravity. 

Building on the seminal works by P. Hofava [1, 2], we have recently proposed a power- 
^ counting renormalizable model for quantum gravity without Lorentz invariance [3] . Remark- 

ably, the model is free of the pathologies [4, 5, 6, 3, 7] present in the original Hofava's proposal 
and associated with the additional mode of the gravitational excitations. This is achieved 
by providing the extra mode with a proper quadratic action around smooth backgrounds. 
We have argued in [3] that this property together with power-counting renormalizability of 
the theory ensures that the theory is weakly coupled all the way up to trans-Planckian ener- 
gies 1 . We also argued in [3] that with appropriate choice of parameters the theory satisfies 



lr To be weakly coupled at all energies the model must fulfill the additional requirement that its marginal 
couplings do not develop Landau poles under the renormalization group flow. In other words, the theory 
must possess a weakly coupled UV fixed point. The Landau poles, if any, appear at exponentially high 
energies and are irrelevant for the purposes of the present paper. 
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phenomenological constraints, and demonstrated this explicitly for the simplest tests pro- 
vided by the gravitational potential between static sources (Newton's law) and homogeneous 
cosmology (Friedmann equation). 

The consistency of the model presented in [3] has been recently questioned in [8], where 
it is claimed that the model suffers from the same kind of strong coupling problem as the 
previous versions of Horava's proposal [6, 7]. The aim of the present note is to show that 
this claim is unfounded. We support our arguments by considering a toy model where the 
absence of strong coupling is demonstrated both using the power counting and at the level 
of scattering amplitudes. We also analyze in more detail the observational bounds on the 
model and determine the window in the parameter space compatible with phenomenology 
and weak coupling. 

Review of extended Hofava gravity 

We start by describing briefly the model [3]. We consider the Arnowitt-Deser-Misner (ADM) 
decomposition for the metric, 

ds 2 = (N 2 - NiN^dt 2 - 2N i dx i dt - -f ij dx i dx j . (1) 

This decomposition defines a foliation of space-time by 3-dimensional space-like surfaces thus 
splitting the coordinates into "space" and "time". We follow [2] and, unlike the case of general 
relativity (GR), consider this foliation structure as physical. This means that the group of 
invariance of the theory is not the full group of 4-dimensional diffeomorphisms, but only its 
subgroup consisting of foliation-preserving transformations 

x H> x(f,x) , t H> t(t) . (2) 

The action of the model is taken in the form 2 

M 2 f 

S = I d 3 xdt^N (KijICi - \K 2 - Vhj, a,]) , (3) 
where Mp is the Planck mass; Kij is the extrinsic curvature tensor 

K i3 = ^ (7« - ViNj - VjNi) , (4) 

2 The 3-dimensional indexes . . are raised and lowered using 7^, and covariant derivatives are associ- 
ated to Jij. 



2 



with the trace K\ 7 is the determinant of the spatial metric 7^; A is a dimensionless constant. 
The "potential" term V[7y,aj] in (3) depends on the 3-dimensional metric and the lapse N. 
The latter enters into the potential through the combination 

at = -jf , (5) 

which is covariant under the symmetry (2). The dependence of the potential term on the 
gradients of the lapse is the key difference of the model (3) compared to the original Hofava's 
proposal [2]. 

The existence of the preferred foliation structure reflects the non-relativistic nature of 
the model: space and time enter into the theory on different footings. This allows to in- 
troduce into the action terms with higher derivatives in spatial directions which improve 
the ultraviolet (UV) behavior of the graviton propagator [2]; at the same time the theory 
remains second order in time derivatives thus avoiding problems with unitarity. Out of the 
previous family of actions (3), one can construct power-counting renormalizable theories by 
considering the scaling transformations [2] 

x i-> & _1 x , t h-> b~ 3 t , (6a) 
N^N , Ni^b 2 Ni , iij h-> iij . (6b) 

Under this scaling, the kinetic part of the action (3) and the operators of dimension 3 6 in 
V are left unchanged (they are marginal) 4 . Operators of lower dimensions in V are relevant 
deformations. According to the standard arguments, considering operators up to dimension 
6 in the potential gives rise to an action which is perturbatively renormalizable. Explicitly, 
the allowed potential term is 

V = — £R — a aid 1 

+ Mp 2 (A 1 AR + A 2 R tj R ij + A^Aa 1 + A^a^) 2 + A^ajR^ + . . .) (7) 
+ M P A (B 1 A 2 R + B 2 R ij R jk R i k + 5 3 a;AV + B^a^f + B h a { a l l aja k R jk + ...), 

where Rij, R are the Ricci tensor and the scalar curvature constructed out of the metric 7^; 
A = 7 lJ 'VjVj, and £, A n , B n are constants. The ellipses represent other possible operators 
of dimension 4 and 6 which can be constructed out of the metric 7^ and are invariant under 



3 We asign dimension —1 to the space coordinate. Then the dimension of time is —3, dimensions of the 
lapse and the 3-dimensional metric are zero, etc. 

4 This is true classically. At the quantum level one expects the coefficients in front of marginal operators 
to acquire logarithmic running under the renormalization group flow. 
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3-dimensional diffeormorphisms 5 . In what follows we set £ = 1 , which can always be achieved 
by a suitable rescaling of time. At low energies the potential is dominated by the operators 
of the lowest dimension, namely, the terms in the first line of (7). This leads to the recovery 
in the infrared of the relativistic scaling dimension —1 for both space and time. 

The explicit breaking of 4-dimensional diffeomorphisms down to the subgroup (2) gives 
rise to the presence of a new scalar gravitational degree of freedom [2, 3, 6]. Its properties 
at the quadratic level were analyzed in [3] where it was shown that the new mode is free of 
pathologies at all energies (it is neither a ghost nor a tachyon) in a wide range of parameters. 
The proper behavior of the mode at low energies is ensured by the following choice of the 
parameters A and a (see Eqs. (3), (7)) 

0<^^, 0<a<2. (8) 

The additional mode does not have a mass gap: at low energies it obeys a linear dispersion 
relation with a velocity generically different from that of gravitons (which is 1 in our choice 
of units). This signals the break down of Lorentz invariance down to arbitrary low energies. 
As we discuss below, this has phenomenological consequences that ultimately set bounds on 
the values of the parameters A and a governing the low-energy physics of the model. 

It is convenient to introduce the covariant form of the theory which we obtain using 
the method described in [6]. One encodes the foliation structure of space-time into a new 
Stiickelberg field <j)(t, x) by identifying the surfaces of the foliation with surfaces of constant 

4>(t, x) = const . (9) 

The invariance of the theory under reparameterization of the foliation surfaces translates 
into the invariance under reparamaterizations of 0, 

<f> H- , (10) 

where / is an arbitrary monotonous function. The quantities appearing in the action (3) 
reduce to the standard geometrical objects (induced metric, extrinsic and intrinsic curvature) 

5 Note that the operators with odd dimensions are forbidden by spatial parity. Similarly, the terms in the 
action with one time derivative of a, are excluded by the time-reversal invariance. We stress that apart from 
these restrictions one must consider all operators of dimension up to 6 and compatible with the symmetries 
(2) to obtain a renormalizable action. Such operators are numerous and only a few of them are written 
explicitly in the above expression. The complete list of terms providing non-equivalent contributions at the 
quadratic level is given in [3] . 
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(12) 



characterizing the embedding of the hypersurfaces defined by (9) in space-time. The central 
object in the construction of these quantities is the unit normal vector 6 u^. Explicitly, 

Note that is automatically invariant under the transformations (10). Other geometrical 
quantities associated to the foliation are constructed out of w M and its derivatives, see [6] for 
details. In this way one obtains the following covariant form of the action (3), 

S = — / d 4 x v /=i| (4) J R+ (A- l)(V/) 2 + a«VV/V„« p 

+ (terms with higher derivatives) | . 

This action describes a scalar-tensor theory of gravity invariant under 4-dimensional dif- 
feomorphisms and the symmetry (10). Furthermore, after fixing the gauge = t, it is 
equivalent to the non-covariant form (3). Thus the covariant (Stiickelberg) formalism makes 
the presence of the extra scalar degree of freedom explicit. 

The first line in (12) contains all the terms with up to two derivatives acting on and 
the metric; it describes the low-energy physics of the model. Note that this low-energy action 
is similar to a special case of the Einstein-aether theory (see [9] for a recent review). The 
difference from the general Einstein-aether theory is that in our case the vector is, by 
its definition (11), hyper surf ace-orthogonal; i.e. it is characterized by a single scalar field 7 . 
Besides the low-energy part, the full action of the model also contains the terms with higher 
derivatives which we do not write explicitly. These terms arise from the second and third 
lines in the potential (7). Their important effect is to modify the dispersion relation of the 
modes at high energies, 



4 6 
22, P , P 



E = cV + ^— + ^— } (13) 
Ml A M$ B 



where E and p are the energy and momentum of the modes, and 

for helicity-2 modes 




(14) 

for scalar graviton. 



6 The Greek indices fi,v,... are raised and lowered using the 4-dimensional metric g^ v , and the covariant 
derivatives with these indices are understood accordingly. 

7 In comparison with [9], we have absorbed one of the free parameters of the most generic action in a 
redefinition of time (£ = 1). 

8 To be precise, the lower expressions holds in the "decoupling limit" when a, |A — 1| <C 1. See [3] for the 
exact expression. 
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For simplicity, we assume c s ~ 1 in what follows. In this case, one reads off the scales 
suppressing the higher derivative terms from (3), (7): 

[a; 1/2 M p , B; l/A M P for helicity-2 modes 

M ** ~ < ' i/2 „/ i/4 ( 15 ) 

I v^^i Mp, oi ' B i Mp for scalar graviton. 

As we now discuss, the presence of these higher-derivative terms is crucial to make the theory 
weakly coupled and renormalizable in the UV. 

Would-be strong coupling and its resolution 

In the covariant language, the issue raised in [8] , can be understood as follows. Let us expand 
the low-energy action of the model (first line in (12)), around the background consisting of 
the Minkowski metric and the Stiickelberg field linearly depending on time, 

9nv(t,x) = V»v + V(^ x ) > 0(*> x ) =< + *(*» x ) • ( 16 ) 
The result has the schematic form 



Ml I (fx 



- hUh - a{da) 2 + (A - 1)(A X ) 2 + (A - l)x(Ax) 2 + • ■ • 



where, for the sake of the argument, we have written down only one of the interaction terms. 
The quadratic part of the perturbed action remains invariant under the relativistic scaling 

x fe _1 x , t b~H , (17a) 
V *-> bh n» > X >->• X ■ (17b) 

The interaction terms for both fields have positive dimensions with respect to this scaling. 
This means that these interactions would become strong at a certain scale A if no new physics 
appeared at a lower scale. Under the assumption (motivated by phenomenological bounds) 
| A — 1 1 ~ ce 1, the covariant formalism with the action (12) allows to readily identify the 
scale A as 



A = y/\X-l\M P ~ ^faMp . (18) 

The scale (18) has been erroneously interpreted in [8] as the UV cutoff of the theory where 
the perturbative description breaks down. Actually, A is only the cutoff of the low- energy 
approximation. In the model described above, the would-be strong coupling is actually 
not present if the higher-derivative operators (which change the scaling dimensions of the 
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interactions) enter into the game at energies lower than (18) (see the related discussion in 
[2]). For this to happen the energy scale of UV physics (15), which we collectively denote by 
M*, must be smaller than A, 9 

M, < v^Mp . (19) 

Then, the power-counting analysis performed in the ADM frame (see (6)) shows that under 
the new scaling the interactions are at most marginal, meaning that there is no strong 
coupling at the scale A. We conclude that the correct interpretation of the scale (18) in the 
model at hand is that of the scale suppressing the higher- derivative operators. 

Let us illustrate our point by a simple toy model. Consider a scalar theory with action 



S = aM 2 P fd 4 xl(ip + J2 a n¥ n ) 

J { n>2 



m A 3 



V , (20) 



where the dimensionless coupling constants a n are assumed to be somewhat smaller than 
1. This scalar theory shares all the relevant properties with our actual gravity theory. At 
low momenta, |A| «C M 2 , the higher derivative terms can be neglected and one obtains the 
following low-energy action 



~>lowE 



aM 2 P I d 4 x| (<p + a nV n ) □</?} • (21) 

J n>2 



Clearly, the invariance of the quadratic part of the action with respect to the relativistic 
scaling transfromations (17a) sets the scaling dimension of tp to be 1. The action contains 
irrelevant interactions under this scaling which naively become strong at the scale 

A = v^Mp . (22) 

However, this is not the case provided M* < A. At momenta above M* the quadratic action 
is dominated by the term with the highest number of spatial derivatives, 



. (23) 

This is invariant under anisotropic scaling transformations (6a) with <p having scaling di- 
mension zero. Consequently, all the interactions in the full action (20) become marginal at 



9 One may worry that the choice of M* (and A) parametrically below Mp introduces a fine-tuning in the 
model. Let us emphasize that this is not the case: having M* well below Mp is technically natural. From 
the point of view of the low-energy theory, the reason is that the cutoff is set by M«, and not Mp. Thus, 
neither Mp nor M* receive large corrections. 
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high energies, and the relative strength of the interaction terms with respect to the free part 
is always small. 

It is instructive to see explicitly how the terms with higher derivatives prevent the theory 
from becoming strongly coupled in the language of scattering amplitudes. A well-known 
manifestation of the breakdown of perturbation theory is the saturation of unitarity bounds 
by tree-level amplitudes (see e.g. [10]). From the low-energy form of the action (20) one 
would conclude that tree-level unitarity is violated at the scale (22) and that perturbation 
theory is no longer valid at higher energies. As we shall now discuss, this conclusion would 
be incorrect. This is essentially due to the peculiar kinematics of theories with anisotropic 
scaling, summarized by the dispersion relation (13), which makes the unitarity bound much 
milder at high energies as compared to the relativistic case. 

To be concrete, we consider the s-channel scattering of two tp quanta with energy E 
in the center of mass frame (which we assume to coincide with the preferred frame). The 
optical theorem yields for the s-wave amplitude 10 M.{2 — > 2), 

2Im.M(2 ->■ 2) = 

j/^^i^^^i'w^-E^O^a:^. (24) 

where the sum runs over all possible final states. Note that this relation holds for arbitrary 
dispersion relation 

E = £{p) . (25) 

Being a sum of positive numbers, the r.h.s of (24) is larger than any of the summands. In 
particular, 

2Im.M(2->2)> J^ 2 ^L- ) \M(2^2)\ 2 S(2E -2£(p)), (26) 

Performing the integrations and using ImM < \M\ one obtains the bound on the absolute 
value of the amplitude, 

\M(2^2)\< 16tt£:'(po)4 > ( 27 ) 

Po 

where po is related to the energy Eq of the incoming particles by the dispersion relation (25). 
This bound simplifies when the dispersion relation is given by a power-law, £{p) = p z /M^ _1 , 

\M\ < 16nz [E /M*f z - 1)/z . (28) 



10 We stick to the 'relativistic' normalization of the 1-particle states |p) = y / 2£(p)a| > |0). This choice leads 
to conventional expressions for the amplitudes at low energies where the relativistic dispersion relation is 
recovered. 



8 



For relativistic particles, z = 1, this takes the familiar form < 16n. However, for 

the case of anisotropic scaling with z > 1 the bound (28) is less restrictive and allows the 
power-law growth of the amplitude with energy. 

Let us check that in the model (20) the bound (27) is satisfied. The dispersion relation 
8{p) = p y/l + (p/M*) 4 interpolates between z — 1 and z = 3 at low and high energies 
respectively. The leading contribution to the tree-level amplitude comes from the diagram 



which is estimated as 



M(2 - 2) ~ , (29) 



where each vertex contributes a factor E^/ y/aMp and the propagator 1/Eq. At low energies, 
p <C M*, the bound (27) reduces to the condition 

p2 <1. 



aM p ~ 

Naively, this would imply the breakdown of tree-level unitarity at ^faMp. However, if 
yfaMp > M*, the low-energy approximation fails and the bound reads instead (for p 3> M*) 

Ml < aM 2 P , (30) 

which is indeed satisfied. Thus, we recover the same result that was derived from the scaling 
analysis (Eq (19)): there are no signals that perturbation theory is breaking down at the 
scale y/aMp. 

A few remarks are in order. The above arguments do not exclude the possibility that 
some marginal coupling of the theory develop a Landau pole when loop corrections are taken 
into account. If this turns out to be the case, the theory will become strongly coupled in the 
deep UV (thus spoiling the UV-completeness of the proposal). Even in this case, though, 
this would happen at an exponentially high energy. For example, in the case of the toy 
model (20) one can estimate this scale as 

^Landau ~ exp 



P (a„(M*))' 

where /3 and 7 are numerical coefficients of order 1. Certainly, the presence or absence 
of Landau poles in the extended Hofava model (3) is an important open issue requiring a 
detailed renormalization group analysis of the theory. 
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Another basic issue concerning the consistency of the theory (3) at the quantum level is 
to demonstrate the absence of anomalies in the symmetry (2). We hope to return to these 
issues in the future. 

Observational bounds on the UV scale 

The weak coupling condition (19) does not allow to take the parameters a, \X — 1| to zero. 
As already emphasized in [3], this implies that the model does not possess a GR limit in 
the IR, since inevitably a gapless scalar polarization persists down to the lowest energies. 
On the other hand, from the remarkable success of GR in the description of low-energy 
gravitational physics, one expects the phenomenological constraints to put upper bounds on 
the parameters a and |A — 1|, and hence on M*. Thus the real physical question is whether 
it is possible to comply with observations without lowering the scale A to an unacceptable 
level. 

From the fact that the low-energy form of the action (12) corresponds to a special case 
of the Einstein-aether theory [9] one expects that the phenomenology of the two models 
may be similar. This expectation is supported by the results [3] for the weak gravitational 
field of static sources (at rest in the preferred frame) and for the expansion of the Universe. 
Similarly to the case of Einstein-aether theory, static sources in the model (3) give rise to a 
linear metric which has the same form as in GR with the Newton's constant 11 

G N = (8ttM|,(1-«/2)) _1 . (31) 

Importantly, this implies that the PPN parameter ^ PPN has its GR value, ^ PPN = 1. The 
cosmological expansion in the model (3) is governed by the standard Friedmann equation with 
the effective gravitational constant G cosm ^ Gn, which again coincides with the situation 
in the Einstein-aether theory. The phenomenological constraint \G cosm /GN — 1| < 0.13 [11] 
sets a mild bound [3]: a, |A — 1| < 0.1. 

Following the guide of the Einstein-aether, one expects the most stringent constraints on 
the model (3) to come from the observational bound on the PPN parameter a PPN which 
characterizes the preferred frame effects due to Lorentz violation (see [12] for the precise 

n The expression (31) is obtained under the assumption that matter couples universally to the metric <? M „. 
A more general situation compatible with low-energy Lorentz invariance in the matter sector is coupling it to 
the universal effective metric g e JJ — + (3u fl u u , where (3 is a dimensionless parameter. This modification 
preserves the GR form of the weak gravitational field but changes the expression for the Newton's constant. 
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definition) 12 . The detailed study of the PPN corrections in the model (3) will be reported 
elsewhere 13 [14]; here we only sketch the estimate for ctf PJV . This parameterizes an angu- 
lar dependent contribution to the Newtonian potential produced by a source moving with 
velocity v with respect to the preferred frame, 

* N = -^(l + ?*!lj sin 2 * 



where m is the mass of the source, and 6 is the angle between the radial vector and the 
velocity of the source with respect to the preferred frame, cos 9 — f • v. From the physical 
point of view, this contribution is due to the interaction via the Lorentz- violating scalar 
mode associated with the vector w M in the action (12) (cf. discussion in [15]). The result 
(31) for static gravitational field shows that for a <C 1 the scalar-exchange amplitude is 
suppressed by a. Thus we conclude that 14 

„PPN „ /QQ\ 

The bound on a P PN following from the observed alignment of the rotation axis of the Sun 
with the ecliptic [12] gives 

a, | A — 1 1 <4x 1(T 7 , 

where we again assume a and |A — 1| to be comparable. This translates into the bound on 
the suppression scale for higher-derivative operators 

M* < 10 15 GeV . (33) 

To our knowledge, this is the strongest upper bound arising from gravitational physics. 

The lower bound on M* from purely gravitational physics is very mild. Direct tests of 
Newton's law at the distances ~ 10 /xm imply [12] 

M*>0.1eV. (34) 

A stronger bound may be obtained under the additional assumption that M* also sets the 
suppression scale for terms with higher powers of momentum p in the dispersion relations of 



12 The constraint due to the absence of gravitational Cherenkov emission by high-energy cosmic rays [13] 
is easily evaded by setting the velocity of the scalar graviton (as well as that of the helicity-2 mode) larger 
or equal than the maximal velocity of matter particles. 

13 The details of the PPN calculations in the model (3) are different from the Einstein-aether case due to 
the absence of the transverse vector mode. 

14 An explicit computation [14] yields af > " PJV = a 2 /2(\ — 1) = a/2c 2 s , where c s is defined in (14). This 
coincides with (32) when c s is of order one. 
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the matter fields, specifically, of photons. Timing of active galactic nuclei [16] and gamma ray 
bursts [17] constrains the value of such terms. Note that odd powers of p can be forbidden 
(at least, in electrodynamics) by imposing parity. Then, the leading contribution to the 
dispersion relation has the form p 4 /M^ which yields [16, 17] 

M* > 10 10 -T- 10 n GeV . (35) 

Let us stress that, unlike the upper bound (18), this lower bound is model dependent: it relies 
on the assumption that the UV modification to the dispersion relation for photons appears 
at the same scale as that for scalar graviton. This need not hold in some formulations of the 
theory. 

It is worth emphasizing the difference between the situation in the model (3) and that in 
Horava's original proposal, both in its projectable and non-projectable versions [2]. In both 
versions, the strong coupling scale calculated within the low-energy effective theory is so low 
that the introduction of any new physics at that scale is phenomenologically unacceptable. 
Indeed, as discussed in [6], in the non-projectable case the strong coupling scale for the 
additional mode is inversely proportional to the curvature radius of the background. It goes 
to zero for flat, cosmological and static backgrounds, invalidating the proposal. 

In the projectable case it was shown [3, 7] that the scalar graviton mode is unstable at 
large wavelengthes. The requirement that the rate of the instability is smaller than the age of 
the Univers (in order not to spoil standard cosmology) gives the bound |A — l) 1 / 2 < H /M*, 
where Hq is the present value of the Hubble parameter and M* is the suppression scale of 
the higher-derivative operators [3, 7]. On the other hand, M* must be smaller than the 
strong coupling scale of the low-energy theory, which in this case is [7] |A — l\Mp. This 
gives the bound M* < (HqMp) 1 ^ 3 ~ (1000 Km) -1 . Comparing this with the experimental 
bound (34), one concludes that the projectable case in the weakly coupled regime is ruled out. 

To sum up, we have shown that the claim [8] about the presence of strong coupling 
problem in the model (3) is unfounded. The absence of strong coupling is actually a built-in 
feature of the model. It suffices for the scale M* suppressing the higher derivative operators 
to be slightly lower than the naive strong coupling scale calculated in the low-energy theory. 
The observational constraints on deviations of gravity from GR place an upper bound M* < 
10 15 GeV. Under the additional (model dependent) assumption that this scale is common for 
gravity and matter sectors one obtains a lower bound M* > 10 10 -r- 10 11 GeV. Within this 
range, to the best of our knowledge, the model is compatible with the existing data. Thus, so 
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far, the model is a phenomenologically viable candidate for a renormalizable quantum theory 
of gravity. Needless to say, whether the theory is truly renormalizable (anomaly free) and 
UV complete remains an important open issue. Another major question is the mechanism 
for the recovery of Lorentz invariance in the matter sector at low energies (see [18, 19] for 
the detailed discussion of the problem). 
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